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Abstract 

Interactions of heavy Major ana neutrinos in the thermal phase of the early universe 
may be the origin of the cosmological matter-antimatter asymmetry. This mech- 
anism of baryogenesis implies stringent constraints on light and heavy Major ana 
neutrino masses. We derive an improved upper bound on the CP asymmetry in 
heavy neutrino decays which, together with the kinetic equations, yields an upper 
bound on all light neutrino masses of 0.1 eV. Lepton number changing processes at 
temperatures above the temperature Tg of baryogenesis can erase other, pre-existing 
contributions to the baryon asymmetry. We find that these washout processes be- 
come very efficient if the effective neutrino mass rh\ is larger than ~ 10~ 3 eV. 
All memory of the initial conditions is then erased. Hence, for neutrino masses 



in the range from J Am s 2 ol ~ 8 x 10~ 3 eV to yj Am^ tm ~ 5 x 10~ 2 eV, which is 
suggested by neutrino oscillations, leptogenesis emerges as the unique source of the 
cosmological matter-antimatter asymmetry. 



1 Introduction 



The explanation of the cosmological baryon asymmetry is a challenge for particle physics 
and cosmology. In an expanding universe, which leads to departures from thermal equilib- 
rium, C, CP and baryon number violating interactions of quarks and leptons can generate 
dynamically a baryon asymmetry [1]. The possible realization of these conditions has first 
been studied in detail in the context of grand unified theories [2, 3]. 

The picture of baryogenesis is significantly changed by the fact that already in the 
standard model of particle physics baryon (B) and lepton (L) number are not conserved 
due to quantum effects [4]. The corresponding non-perturbative AB = 3 and AL = 3 
processes are strongly suppressed at zero temperature. However, at temperatures above 
the critical temperature Tew of the electroweak transition they are in thermal equilibrium 
[5] and only the difference B — L is effectively conserved. 

During the past years data on atmospheric and solar neutrinos have provided strong 
evidence for neutrino masses and mixings. In the seesaw mechanism [6] the smallness 
of these neutrino masses m v is explained by the mixing mc of the left-handed neutrinos 
with heavy Majorana neutrinos of mass M, which yields the light neutrino mass matrix 

m v = -m D ^m T D . (1) 

Since m D = 0(v), where v ~ 174 GeV is the electroweak scale, and M >> v, the neutrino 
masses m v are suppressed compared to quark and charged lepton masses. CP violating 
interactions of the heavy Majorana neutrinos can give rise to a lepton asymmetry and, 
via the AB = 3 and AL = 3 sphaleron processes, to a related baryon asymmetry. This 
is the simple and elegant leptogenesis mechanism [7]. 

Leptogenesis is a non-equilibrium process which takes place at temperatures T ~ M ± . 
For a decay width small compared to the Hubble parameter, Ti(T) < H(T), heavy 
neutrinos are out of thermal equilibrium, otherwise they are in thermal equilibrium. A 
rough estimate of the borderline between the two regimes is given by r x = H(M 1 ) (cf. [8]). 
This is equivalent to the condition that the effective neutrino mass rh\ = (m\ } mD)ii/Mi 
equals the 'equilibrium neutrino mass' 



16tt 5 / 2 , /0 V 



n; gl/*-£-~10- 3 eV , (2) 



3^ * M 1 



P i 



where we have used M p i = 1.2 x 10 19 GeV and g* = 434/4 as effective number of degrees 
of freedom. For rh\ > m* (mi < m*) the heavy neutrinos of type Ni are in (out of) 
thermal equilibrium at T = M x . 
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It is very remarkable that the equilibrium neutrino mass m* is close to the neutrino 
masses suggested by neutrino oscillations, a/ Am s 2 ol ~ 8 x 1CT 3 eV and a/ Am^ tm ~5x 
1CT 2 eV. This suggests that it may be possible to understand the cosmological baryon 
asymmetry via leptogenesis as a process close to thermal equilibrium. Ideally, AL = 1 
and AL = 2 processes would be strong enough at temperatures above M 1 to keep the 
heavy neutrinos in thermal equilibrium and weak enough to allow the generation of an 
asymmetry at temperatures below Mi. 

An analysis of solutions of the Boltzmann equations shows that this is indeed the 
case if light and heavy neutrino masses lie in an appropriate mass range. In general, 
the final baryon asymmetry is the result of a competition between production processes 
and washout processes which tend to erase any generated asymmetry. Unless the heavy 
Majorana neutrinos are partially degenerate, M2$ — Mi < Mi, the dominant processes are 
decays and inverse decays of Ni and the usual off-shell AL = 1 and AL = 2 scatterings. 
The final baryon asymmetry then depends on just four parameters [9] : the mass Mi of 
Ni, the CP asymmetry Si in N\ decays, the effective neutrino mass fhi and, finally, the 
sum of all neutrino masses squared, m 2 = ml + m\ + m\, which controls an important 
class of washout processes. Together with the two mass squared differences Am 2 atra and 
Ara^, the sum m 2 determines all neutrino masses. Using an upper bound on the CP 
asymmetry Si [10, 11], an upper bound on all light neutrino masses of 0.2 eV has recently 
been derived [12]. 

In this paper we extend the previous analysis in two directions. We derive an improved 
upper bound on the CP asymmetry which leads to a more stringent upper bound on light 
neutrino masses. In addition, we study in detail the washout of a pre-existing B — L 
asymmetry, which yields a lower bound on the effective neutrino mass rhi. In this way 
we obtain a window of neutrino masses for which leptogenesis can explain the observed 
cosmological baryon asymmetry, independent of initial conditions. 

The paper is organized as follows. In Section 2 we derive an improved upper bound 
on the CP asymmetry e\ and illustrate how it can be saturated for specific neutrino 
mass matrices. Theoretical expectations for the range of neutrino masses are discussed 
in Section 3. In Section 4 we then derive upper bounds on the light neutrino masses in 
the cases of normal and inverted hierarchy, and we discuss the stability of these bounds. 
Section 5 deals with the washout of a large initial B — L asymmetry, and a summary of 
our results is given in Section 6. 
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2 Bounds on the CP asymmetry 



Given the masses of heavy and light Majorana neutrinos the CP asymmetry S\ in the 
decays of Ni, the lightest of the heavy neutrinos, satisfies an upper bound [10, 11]. In 
the following we shall study under which conditions this upper bound is saturated and 
how it depends on the effective neutrino mass rh\ which plays an important role in the 
thermodynamic process of leptogenesis. 

The standard model with right-handed neutrinos is described by the lagrangian, 

C m = hijluVRjcfr + T^MijP^VRj + h.c. , (3) 

where M is the Majorana mass matrix of the right-handed neutrinos, and the Yukawa 
couplings h yield the Dirac neutrino mass matrix mo = hv after spontaneous symmetry 
breaking, v = ((f)). We work in the mass eigenstate basis of the right-handed neutrinos 
where M is diagonal with real and positive eigenvalues Mi < M 2 < M 3 . The seesaw 
mechanism [6] then yields the light neutrino mass matrix 

m u = -m D ^m T D , (4) 

which can be diagonalized by a unitary matrix 



/mi 

m 2 | . (5) 
\ 




with real and positive eigenvalues satisfying m\ < m 2 < m 3 . 

It is convenient to work in a basis where also the light neutrino mass matrix is diagonal. 
In this basis the Yukawa couplings are 

h = U^h . (6) 

As a consequence of the seesaw formula the matrix fl, 

v 

^ij = / fj hij , (7) 

\JmiMj 

is orthogonal, f2f2 T = f2 T fi = / [13]. It is then easy to show that the CP asymmetry e\ 
[14]-[16] is given by (cf., e.g., [9]) 



where Am^ = mf — m\. 

The CP asymmetry e 1 is bounded by the maximal asymmetry e™ ax [12], 

I I < r max _ 3 Mi (Aml tm + Am 2 sol ) 

l£ll ~ £l "16^ m 3 • (9) 

As we will now show, this bound holds for arbitrary values of m 2 , i.e. for normal and for 
inverted hierarchy, and it is saturated in the limit mi — > 0. 
Consider the normalized Yukawa couplings 

h 2 

' : ■>■, I ill, ■ (10) 



with 

o<N<i, ^2\zi\ = i. (ii) 

i 

The orthogonality condition (f2 T fi)n = 1 yields the additional constraint 



In the new variables the CP asymmetry reads 



3 M 1 /Ami! , Am|i , 
£ i = 77 T V2 + -2/3 • (13) 



167T v 2 \ m 2 m 3 

Since m 3 > m 2 , one also has Am|i/m 3 > Am 21 /m 2 . This suggests that the maximal CP 
asymmetry is reached for maximal y 3 . 

Suppose now that 1 — y 3 = 0(e). Because of Eqs. (11) this implies 1/2, yi and all X{ 
have to vanish in the limit e — > 0. The orthogonality condition (f2 T f2)n = 1 yields 

^ + ^ + ^=0, (14) 
mi m 2 m 3 

— xi H rr 2 H rc 3 = 1 . (15) 

mi m 2 m 3 

Since m 2 > 0, these conditions are satisfied for maximal y 3 , if y 2 = x 2 = x 3 = and 

mi , y 1 oc e , (16) 

ii oc e a , xi oc e 1_a , < a < 1 . (17) 

Note that in the limit e — > 0, iVi couples only to / 3 0. For a > 0, iVi decouples completely, 
since /i 2 x = {h)h)nZi and {h)h)u oc m-i. 



An explicit example, which illustrates this saturation of the CP bound, is given by the 
following orthogonal matrix, 



( A 





1 




-B \ 


A ) 



with 



B 



-be a , A 



m 3 M x 

The corresponding Yukawa couplings squared are 



B' 



b>0. 



One obviously has x 2 = x 3 
matrix of Yukawa couplings, 

h = 



V2 



m 3 

0, and y 3 — > 1, Xi,yi — > in the limit e 



(18) 

(19) 

(20) 
0. The 



v A mo, 





V 



-Vibe" 





v / m-2-M2 

V 





„2 



l M l Ut 



(21) 



becomes diagonal in the limit e — > for a > 0. Hence, in this basis, the large neutrino 
mixings are due to the charged lepton mass matrix. 

This example illustrates that rh\ can be arbitrary in the limit m\ — > 0. It approaches 
b 2 v 2 /Mi for a = 0, while it goes to for a > 0. Hence, the maximal CP asymmetry 
(9) can be reached for arbitrary values of ni2 and rh\. For a given CP asymmetry, the 
maximal baryon asymmetry is reached in the limit rh\ — > 0, assuming thermal initial Ni 
abundance. The corresponding, model independent lower bound on the heavy neutrino 
mass Mi was determined in [9] to be > 4 x 10 s GeV. If the Yukawa couplings h are 
restricted, a more stringent lower bound on Mi can be derived [17]. 

The above discussion can easily be extended to derive the maximal CP asymmetry in 
the case of arbitrary rhi. Since m 3 > m 2 > mi, one again has x 3 = x 2 = t/2 = 0. From 
Eqs. (14), (15) one then concludes 

mi mi 

yi = 1/3 , xi = — . 

m 3 mi 



(22) 



Together with the constraint (cf. (11)), \fx\ + y\ + |y 3 | = 1, these conditions determine 
1 2/3 1 as function of mi, m 3 and fh 1 . Inserting the result into Eq. (13) yields the improved 
upper bound 



3 M x m 3 

167T V 2 



1 - 



mi 
m 3 



m% 



1/2' 



(23) 
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For mi = the result coincides with the previous bound (9). For < m x < rh\ the 
new bound is more stringent. In particular, e™ ax = for rh\ — m\. Note that according 
to Eq. (23) the only model independent restriction on the effective neutrino mass is 
rhi > m\. The improved upper bound on the CP asymmetry implies also a bound on the 
light neutrino masses which is more stringent than the one obtained in [12]. This will be 
discussed in Section 4. 

3 Range of neutrino masses 

At present we know two mass squared differences for the light neutrinos, which are deduced 
from the measurements of solar and atmospheric neutrino fluxes. In addition we have 
information about elements of the mixing matrix U in the leptonic charged current. Since 
U could be entirely due to mixings of the charged leptons, this does not constrain the light 
neutrino mass matrix in a model independent way. The light neutrino masses mi < m 2 < 
m3 can be either quasi-degenerate or hierarchical, with m 2 — vn\ <C — m 2 ('normal 
hierarchy') or m 3 — m 2 <C m 2 — mi ('inverted hierarchy'). The best information on the 
absolute neutrino mass scale comes from neutrinoless double /3-decay, which yields an 
upper bound on the light Majorana neutrino masses of about 1 eV [18, 19]. 

A crucial quantity for thermal leptogenesis is the effective neutrino mass rhi which is 
always larger than mi [20], as one easily sees from the orthogonality of Q (cf. (7)), 

v 2 

i i 

> mi l^al > ^i5Z Re ^i) = mi • ( 24 ) 

i i 

As we saw in the previous Section, the maximal CP asymmetry is reached for mi = 0, 
such that m 2 ~ a/ Am^ and m 3 ~ y/ Am^. 

There is no model independent upper bound on rh\. However, if there are no strong 
cancelations due to phase relations between different matrix elements, one has 

< ra>3 ^ l^ql ^ ^1 ^ — 7^3 ■ (25) 

i i 

Hence, the natural range for the effective neutrino mass is m x < rh\ < m 3 . In fact, we 
are not aware of any neutrino mass model where this is not the case. 

It is instructive to examine the range of rhi also in the special case = e™ ax . As we 
saw in the previous section this case is realized for y 2 = x 2 = x 3 = 0, corresponding to 
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Re(0| 1 ) = Re(Jl| 1 ) = Im^U = 0. The orthogonality condition then implies In^Q^) = 
— Im^g-J and Re(r2^) = 1. Hence, for maximal CP asymmetry one has 



showing that the value of mi is tuned by just one quantity. For Im^U = 0, one 
has mi = mi, while the case rh\ ^> m 3 corresponds to a fine tuned situation in which 



If the observed large mixing angles in the leptonic charged current originate from the 
neutrino mass matrix, which appears natural since their Majorana nature distinguishes 
neutrinos from quarks, the masses m\ and rh\ are related to m 2 and 777.3. The seesaw 
mechanism together with leptogenesis then also constrains the heavy Majorana neutrino 
masses. 

Large mixing angles are naturally explained if neutrino masses are quasi-degenerate 
[23]. One then has rfi2 ~ ITI3 > 0.1 eV. However, as shown in [9, 12] and 

further strengthened in the following Section, quasi-degenerate neutrinos are strongly 
disfavored by thermal leptogenesis. A possible exception is the case where also the heavy 
Majorana neutrinos are partially degenerate. One then gets an enhancement of the CP 
asymmetry which allows one to increase the neutrino masses and still have successful 
leptogenesis. Models with AM 2 i/M 1 = (M 2 - M 1 )/M 1 < 5 x 10~ 2 and AM 2] A/, = 
5 x 1CT 7 have been considered in refs. [25] and [26], respectively. Note, however, that in 
these examples the light neutrino masses are not quasi-degenerate. We shall pursue this 
case further in Section 4.3. 

The neutrino mass pattern with inverted hierarchy has also received much attention in 
the literature. There is, however, the well known difficulty of this scenario to fit the large 
angle MSW solution [27, 28]. We also do not know any model with inverted hierarchy 
which incorporates successfully leptogenesis, and we shall therefore not pursue this case 
further. 

We are then left with the case of neutrino masses with normal hierarchy. There are 
many neutrino mass models of this type with successful leptogenesis. The mass hierarchy 
is usually controlled by a parameter e C 1. For the effective neutrino mass one can then 
have, for instance, mi ~ m 2 (cf. [26, 29]). A simple and attractive form of the light 
neutrino mass matrix, which can account for all data, is given by [30, 31], 




(26) 



\lm(Q 2 31 )\ = l Im (^n)l > Re(^ii) = 1 - 



( 



e 



e 



\ 




e 1 1 



(27) 



Veil/ 
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where coefficients 0(1) have been omitted. This form could follow from a U(l) fam- 
ily symmetry [32] or a relation between the hierarchies of Dirac and Majorana neutrino 
masses [33]. In the second case one has mi, mi ~ em 3 and rh\ ~ m%, which is compatible 
with leptogenesis. The structure of the mass matrix (27) as well as predictions for the 
coefficients 0(1) can be obtained in seesaw models where the exchange of two heavy Ma- 
jorana neutrinos dominates [34]. In all these examples the range of the effective neutrino 
mass is mi < fhi < m 3 . 

Thermal leptogenesis also leads to a lower bound on Mi, the smallest of the heavy 
neutrino masses [35, 11]. In the minimal scenario, where the heavy neutrinos are not 
degenerate, one obtains the lower bound Mi > 4 x 10 8 GeV [9]. It is reached for maximal 
CP asymmetry (mi = 0), minimal washout (mi — > 0), and assuming thermal initial N\ 
abundance. The bound becomes more stringent for restricted patterns of mass matrices 
[17]. It can be relaxed if the heavy neutrinos are partially degenerate [24, 25, 26]. 



4 Improved upper bounds on neutrino masses 

4.1 Maximal asymmetry and CMB constraint 

It is useful to recast the maximal CP asymmetry (23) in the following way, 

e™(Mi, mi ,m) = 10- (J^) ^ ftm^m) , (28) 



where m atm = y/Aml tm + Am 2 sol , m = (16tt/3) (v 2 /10 10 GeV) ~ 0.051 eV, and 

(m 3 - mi J 1 + 

p(mi,m)= ± y - ^-L<\ . (29) 

m atm 

The maximal value, f3 — 1, is obtained for mi = 0. Note, that m atm = m for the 
best fit values extracted from the KamLAND data [22], Am^, = 6.9 x 10~ 5 eV 2 , and the 
Super Kamiokande data [21], Am 2 tm = 2.5 x 10~ 3 eV 2 . 

We will calculate particle numbers and asymmetries normalized to the number of 
photons per comoving volume before the onset of leptogenesis at [9]. For zero initial 
B — L asymmetry, i.e. N l B _ L = 0, the final B — L asymmetry produced by leptogenesis 
is given by 

^-L = -|ei«f> (30) 

where Kf is the 'efficiency factor'. In the minimal version of thermal leptogenesis one 
considers initial temperatures T; > Mi, where M 1 is the mass of the lightest heavy 
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neutrino N±. In this case Kf < 1, and the maximal value, K{ — 1, is obtained for thermal 
initial iV\ abundance in the limit rh\ — > 0. The heavy neutrinos JVi then decay fully out of 
equilibrium at temperatures well below Mi, producing a. B — L asymmetry which survives 
until today since all washout processes are frozen at temperatures T <C M\. 

In the case of general initial conditions and arbitrary values of fhi, the efficiency factor 
Kf has to be calculated by solving the Boltzmann equations [36, 37, 38, 39, 9], 

^ = -(D + S) (N Nl - N%) , (31) 

= - ei D(N Nl -N^)-WN B . Ll (32) 

where z = M\jT . There are four classes of processes which contribute to the different 
terms in the equations: decays, inverse decays, AL = 1 scatterings and processes mediated 
by heavy neutrinos. The first three all modify the N\ abundance. Denoting by H the 
Hubble parameter, D = Yd/ {Hz) accounts for decays and inverse decays, while S = 
Yg/{H z) represents the AL = 1 scatterings. The decays are also the source term for 
the generation of the B — L asymmetry, the first term in Eq. (32), while all the other 
processes contribute to the total washout term W = Y W /{H z) which competes with the 
decay source term. 

We take into account only decays of Aq, neglecting the decays of the heavier neutrinos 
N 2 and N 3 . These decays produce a B — L asymmetry at temperatures higher than Mi. 
As we shall see in Section 5, the washout processes at T ~ Mi very efficiently erase 
any previously generated asymmetry. Even in the case of very small mass differences the 
decays of N 2 and ./V3 do not change significantly the bound on the light neutrino masses, 
which is our main interest. This will be discussed in Section 4.3. 

f 

B—L 



The baryon to photon number ratio at recombination, tjb, is simply related to iV 



by r\B = { a l f) N { B _ L , where a = 28/79 [40] is the fraction of B — L asymmetry which is 
converted into a baryon asymmetry by sphaleron processes, and / = N™ c /N* = 2387/86 
accounts for the dilution of the asymmetry due to standard photon production from the 
onset of leptogenesis till recombination. ^ ax , the final baryon asymmetry produced by 
leptogenesis with maximal CP asymmetry, i.e. Si = e™ x ), is given by 

?7 max ~ 0.96 x 10~ 2 er x ^f • (33) 

This quantity has to be compared with measurements of the CMB experiments 
BOOMerANG [41] and DASI [42], 

V C B MB = (6.01S;J) x 10- 10 . (34) 
10 



The CMB constraint then requires 7]^ x > rf^ , and we will adopt for ry^ the 3a 
lower limit, {t]^ mb ) Xow = 3.6 x 1(T 10 . 

In [9] we showed that rf^, 3 * depends just on the three parameters mi, Mi and m. 
Thus, for a given value of m, the CMB constraint determines an allowed region in the 
(fhi, Mi)-plane. It was also shown that there is an upper bound for m above which no 
allowed region exists. In [12], based on the bound (9) for the CP asymmetry, m < 0.30 eV 
was derived as upper bound on the neutrino mass scale. In the following we shall study 
the allowed regions in the (fhi, Mi)-plane for different parameters m using the improved 
bound on the CP asymmetry (23) and in this way determine a new improved bound on 
m. 



4.2 Numerical results 

The neutrino masses mi and m 3 depend in a different way on m in the cases of normal 
and inverted hierarchy, respectively. Hence, also the dependence of the function (3 on m is 
different for these two mass patterns. This leads to different maximal baryon asymmetries 
77g ax , and therefore to different upper bounds on m, in the two cases which we now study 
in turn. 

For neutrino masses with normal hierarchy one has 

m l ~ m l = A ™atm > m 2 - m l = Am sol , (35) 

and the dependence on m is given by 

ml = - (m 2 + 2Am 2 tm + Am 2 ol ) , (36) 

™l = ^K-AmL + Ami,) , (37) 

ml = \{rn 2 - Am 2 tm - 2Am s 2 ol ) . (38) 

These relations are plotted in Fig. 1. Note that there is a minimal value of m, correspond- 
ing to mi = 0, which is given by m min = a/ Am 2 tm + 2 Am 2 ol ~ 0.051 eV. 

Fig. 2 shows the lines of constant maximal baryon asymmetry rf^** = (r/ ( ^ MB )i ow (thick 
lines) and rfg** = 10 -10 (thin lines) in the (fhi, Mi) -plane for different choices of m and 
assuming zero initial Ni abundance. The allowed regions (the filled ones) correspond to 
the constraint > (r]B MB )\ ow . The largest allowed region is obtained for m = m min , 
since in this case the CP asymmetry is maximal, i.e. (5 — 1 for any value of fhi, and the 
washout is minimal. Note that a different choice for the initial N x abundance would have 
affected the final baryon asymmetry only for fhi < m*. The case of an initial thermal 

11 




m min 0,1 0.20 eV 

rfi/eV 



Figure 1: Neutrino masses as functions of m for normal hierarchy (cf. (36)- (38)). 
10 5 10 4 10 3 10 2 10 1 10° 




m/eV 

Figure 2: Normal hierarchy case. Curves, in the (mi-Mi)-plane, of constant ??^q x = 10 -10 
(thin lines) and tj^q x = 3.6 x 1CP 10 (thick lines) for the indicated values of m. The filled 
regions for tj^q x > 3.6 x 10 -10 are the allowed regions from CMB. There is no allowed region for 
m = 0.20eV. 
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abundance has been studied in [9] . When m increases different effects combine to shrink 
the allowed region until it completely disappears at some value m max . 

We have determined this value with a numerical uncertainly of 0.01 eV. From Fig. 2 
one can see that there is a small allowed region for m = 0.19 eV, whereas we found no 
allowed region for m = 0.20 eV. Hence, the value of m max is somewhere in between and 
we can conclude that in the case of normal hierarchy, 

m<0.20eV. (39) 

Using the relations (36)-(38), one can easily translate this bound into upper limits on the 
individual neutrino masses (cf. Fig. 1), 

mi,m 2 < 0.11 eV , m 3 <0.12eV. (40) 

The case of an inverted hierarchy of neutrino masses corresponds to 

m 2 3 - m 2 2 = Am 2 sol , m 2 2 - m\ = Am 2 tm , (41) 

and the relations between the neutrino masses and m are 

ml = \ (m 2 + Am 2 tm + 2 Am 2 ol ) , (42) 

ml = \ (m 2 + Am 2 tm - Am 2 ol ) , (43) 

ml = \{m 2 -2 Am 2 tm - Am s 2 ol ) . (44) 

We have plotted these relations in Fig. 3. The minimal value of m, corresponding to 
m 1 = 0, is now m min = a/ 2 Am 2 tm + Am 2 ol ~ 0.072 eV. 

The curves of constant rf^ &yi are shown in Fig. 4 for different values of m. The largest 
allowed region is again obtained for m = m min . One can see that this time there is a tiny 
allowed region for m = 0.20 eV and no allowed region for m = 0.21 eV. Therefore, in the 
case of inverted hierarchy the upper bound is slightly relaxed, 

m<0.21eV. (45) 

Using the relations (42)-(44) one can again translate the bound on m into bounds on the 
individual neutrino masses, 

m 1 < 0.11 eV , m 2 , m 3 < 0.12 eV . (46) 

Let us now discuss the different effects which combine to shrink the allowed region 
when the absolute neutrino mass scale m increases, thus yielding the upper bound. The 
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Figure 3: Neutrino masses as functions of m for inverted hierarchy (cf. (42)- (44)). 
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Figure 4: Inverted hierarchy case. Curves, in the (mi-Mi)-plane, of constant f?^Q X = 10 -10 
(thin lines) and rf^^ = 3.6 x 1CP 10 (thick lines) for the indicated values of m. The filled 
regions for tj^q x > 3.6 x 10 -10 are the allowed regions from CMB. There is no allowed region for 
m = 0.20eV. 
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Figure 5: The CP asymmetry global suppression factor /3 max for normal and inverted hierarchy. 

first effect is that away from the hierarchical neutrino case, for to > m m i n and mi > 0, the 
maximal CP asymmetry reduces considerably. This can be seen in terms of the function 
(3 (cf. (29)) which is conveniently expressed in the form 

P = Pm^(m)f(m 1 ,m) . (47) 

The first factor, /3 max = (m 3 — mi)/m atm = m at m/(^3 + wii), is the maximal value of 
(3 for fixed to; /3 max decreases oc 1/m for m ^> m min (cf. Fig. 5). This implies that for 
increasing m there is an overall suppression of the maximal baryon asymmetry in the 
whole (nil, Mi)-plane [11]. In particular the lower limit on M\ becomes more stringent. 

The factor /(mi, to) = 1, for any value of mi, if to = m min (mi = 0). In the case 
m > m min (mi > 0) it vanishes for rh\ = m 1 and grows monotonically to 1 with increasing 
TOi (cf. Fig. 6). Thus for m x > the function / gives a further suppression of the CP 
asymmetry, in addition to the one from /9 max < 1. This suppression is strong for rh\ > TOi 
and disappears for to-i ^> mi. Hence the decrease of the maximal CP asymmetry for 
m > m min shrinks the allowed region most at small mi > mi and at small M\. Note that 
the difference between the allowed regions for normal and inverted hierarchy is accounted 
for by the different values of (5 for a given value of m. In the case of inverted hierarchy 
(3 is larger for any value of fh\ and m > m™ m (cf. Figs. 5,6). The effect is maximal for 
m = m™[ n where (3 mv = 1 while (3 nor ~ 0.6. For larger values of m > m™^ n , and also 
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Figure 6: The function /(mi, to = 0.15 eV) for normal (inverted) hierarchy. It is defined for 
toi > mi ~ 0.08 (0.07) eV. 

fhi ^> m I 1 lor , the ratio / 3 mv / / 9 nor becomes very close to 1. This situation is realized when 
m approaches its upper bound. This explains why the upper bound on m is only slightly 
relaxed in the case of inverted hierarchy. 

The second effect, which shrinks the allowed region when m increases, is the enhance- 
ment of washout processes. In [9] we showed how the total washout rate can be written 
as the sum of two terms, (W — AW) oc rh\ and AW oc M 1 m 2 . The first term is respon- 
sible for the reduction of the allowed region at large rh\. The second term leads to the 
boundary at large M\ . The combined effect shrinks the allowed region with increasing m 
at large Mi and at large rh\. 

One can see how this second effect reduces the allowed region, independent of the 
maximal CP asymmetry decrease, by comparing the two largest allowed regions for normal 
and inverted hierarchy; they correspond to the two different values of m min (cf. Fig. 2 
and Fig. 4). Since f3 — 1 in both cases, the entire difference is due to the different 
washout effects. They are larger in the case of inverted hierarchy because m min is about 
~ y/2 higher than in the normal hierarchy case. One can see how, for a fixed value of 
fhi, the maximal value of M 1 is approximately halved in the inverted hierarchy case. 
Correspondingly, the maximal allowed value of fhi is lower for inverted hierarchy than for 
normal hierarchy. 
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In summary, within the theoretical uncertainties, leptogenesis cannot distinguish be- 
tween normal and inverted hierarchical neutrino mass patterns. However, our new analysis 
confirms and strengthens the results of [9, 12] that quasi-degenerate neutrino masses are 
strongly disfavored by leptogenesis, by putting the stringent upper bound of 0.12 eV on 
all neutrino masses. 



4.3 Stability of the bound 

The numerical results can be very well reproduced analytically [43] . This procedure is not 
only able to yield the correct value of m max but also reveals some general features which 
in the numerical analysis may appear accidental. 

For m = m max , at the peak value of maximal asymmetry, such that ?? max = r/g MB , one 

has 

( \ 

mi La* = m max + O -2^ , (48) 



™max 



M!| max ~ 1.6 x 10 13 GeV ( ) . (49) 



(L2eV\ 



2 



777 



max 



The value of m max is slightly different for normal and inverted hierarchy, respectively, 

«°ax) 2 = (<ax) 2 - ^Ln + ^K m Kj , (50) 

«Tax) 2 = «ax) 2 + lml tm + OimiJK.J , (51) 
where m max is the zero-th order approximation. This implies 

KL) 2 - «Tx) 2 = mi tm + 0{mi t Jr?tJ • (52) 

Besides gaining more insight into the numerical results, the analytic procedure also 
allows to find the dependence of the bound on the involved physical parameters and to 
study in this way its stability. 

Consider first the dependence on the experimental quantities rf^ vlB ', Am^ tm and Am^. 
Since Am^ <C Am atm , the dependence on Am s 2 ol is so small that it can be neglected, 
yielding m atm ~ a/ Am atm . The analytic procedure shows that m max oc (fiatm/'JB™) 1 ^- 
From the numerical result, found for rf^ MB = 3.6 x 1CT 10 and m atm = m ~ 0.051 eV, and 
one then obtains in general 



i ,i 



- max - 0.175 eV [—) • (53) 
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Using Eq. (53) one immediately gets the central value of m max . Note also that for 
^cmb _ iQ-w anc j matm — mQ; one m^ nax ~ 0.275 eV. This is confirmed by the 
numerical results. We still find iso-lines i]^ MB = 1CT 10 for m = 0.27 eV, whereas this is 
not the case anymore for m = 0.28 eV. From Eq. (53) one obtains as estimate for the 
relative error, 

Sm max = \ (S V C B MB + 6ml tm ) . (54) 

According to Eq. (34) the la standard error on r/ B MB is about 15% while om 2 &tm ~ 25% 
[21]. We thus obtain 5m max ~ 10%, which corresponds to the absolute error Am max ~ 
0.02 eV. In the coming years the errors on rf^ MB and m atm will be greatly reduced by 
the satellite experiments MAP [44] and Planck [45], and by the long baseline experiments 
Minos [46] and CNGS [47], respectively, and consequently the error on m max will be 
considerably reduced. 

Another important question concerns the enhancement of the maximal CP asymmetry 
when AM 2 i = M 2 — M l5 where M x and M 2 are the masses of the heavy neutrinos Aq and 
N 2 , becomes comparable to or smaller than M x itself. As long as the mass splitting is 
larger than the decay widths, the enhancement is given by [15, 16], 



2 

= - x 



[1 + x) In 



1 + x\ 2 — x 



x J 1 — x 



(55) 



where x = (M 2 /Mi) 2 . Note, that £ approaches 1 for x ^> 1. The value of m max increases 
like £ 1//4 [43], and it is therefore easy to see how the bound on m gets relaxed for small 
values of the mass difference AM 2 i- 

In Fig. 7 we have plotted the enhancement £ — 1 and the central value of m max , together 
with its la limits, as function of AM 21 /M X . For AM 21 /M 1 > 1 the bounds (39), (45) are 
recovered. Only for values AM 2 i/M 1 < 0.1 the bound gets relaxed in an appreciable 
way. An increase of m max by a factor ~ 3, allowing quasi-degenerate neutrino masses of 
0.4 eV, which could be detected with the KATRIN experiment [48], requires degeneracies 
AM 2 i/Mi, AM 31 /Mi < 10~ 3 . 

In the regime AM 2 i < Mi also decays of N 2 have to be taken into account. As we 
shall see in the next section, for larger mass splittings an asymmetry generated in N 2 
decays would be washed out before T ~ Mi, and it is then sufficient to consider only Aq 
decays. However, even for AM 21 /M 1 < 0.1, it is easy to see that the effect of such an 
additional asymmetry on the bound is small compared to the effect of the CP asymmetry 
enhancement described above. The largest effect would be obtained for s™ 3 * = £ max and 
fh 2 <mi, corresponding to a doubled heavy neutrino abundance without any washout 
enhancement. In this extreme case the bound is relaxed at most by a factor 2 1 / 4 ~ 1.2. 
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Figure 7: The CP asymmetry enhancement £ — 1 (short dashed line) and (m^^} ± Am^ iax for 
normal hierarchy (solid and dashed lines) as functions AM21/M1. 

Even for three degenerate neutrinos, with both AM21/M1 <C 1 and AM 31 /M! <C 1, 
the effect could relax the bound not more than by a factor 3 1 / 4 ~ 1.3. Hence, the CP 
enhancement represents the dominant effect and we can conclude that the bound on m can 
only be evaded in case of an extreme degeneracy among the heavy Majorana neutrinos. 

A further important issue is the effect of supersymmetry on the bound. In this case 
the maximal CP asymmetry is about twice as large which could relax the bound by a 
factor 2 1//4 ~ 1.2. However, washout processes are also considerably enhanced [38]. This 
effect goes into the opposite direction and is actually stronger, so that one can expect a 
slightly more stringent bound on to. A detailed calculation will be presented in [43]. 

We conclude that the leptogenesis upper bound on neutrino masses is very stable. 
The essential reason is that, at to = m max , the peak value rf^ 3 * oc 1/m 4 ^. Hence, any 
variation of the final asymmetry results into change of m max which is almost one order 
of magnitude smaller. The same argument applies also to the theoretical uncertainties. 
Although the various corrections to the Boltzmann equations still remain to be calculated, 
we do not expect a relaxation of TO max by more than 20%. In fact, we expect that the 
corrections will go in the direction of lowering the prediction on the final asymmetry, 
which will make the bound on to more stringent. 

For particular patterns of neutrino mass matrices stronger bounds on the light neutrino 
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masses can be obtained. For instance, one can study how the upper bound changes if 
Mi is required to be smaller than some cut-off value M*. For M{ > Mi| max ~ 10 13 GeV 
(cf. (49)) the bound does not change. For smaller values of Mi the bound becomes more 
stringent. For example, from Figs. 2 and 4 one can see that the cut-off Mi < 5 x 10 12 GeV 
leads to the bound m < 0.15 eV, which corresponds to mi < 0.08 eV. For a restricted 
mass pattern, and neglecting AW oc Mi m 2 washout terms, less stringent bounds have 
been found in [17] for the same cut-off value of Mi. 



5 Dependence on initial conditions 

A very important question for leptogenesis, and baryogenesis in general, is the dependence 
on initial conditions. This includes the dependence on the initial abundance of heavy 
Majorana neutrinos, which has been studied in detail in [9], and also the effect of an 
initial asymmetry which may have been generated by some other mechanism. In the 
following we shall study the efficiency of the washout of a large initial asymmetry by 
heavy Majorana neutrinos. 

For simplicity, we neglect the small asymmetry generated through the CP violating 
interactions of the heavy neutrinos, i.e. we set Si = 0. The kinetic equation (32) for the 
asymmetry then becomes 

^^ = -WN B . L , (56) 

where —N B _ L is the number of lepton doublets per comoving volume. The final B — L 
asymmetry is then given by 

N B _ L = u( Zi )N B _ L , (57) 

with the washout factor 

uj( Zi ) = e-^ dzW{z) . (58) 

In Eq. (56) W(z) = Fw(z)/H(z)z is the rescaled washout rate, where H(z) is the 
temperature-dependent Hubble parameter, receives contributions from inverse decays 
(F ID ), AL = 1 processes (r 0)t , F^ s ) and AL = 2 processes (T N , F N>t ) (cf. [9]), 

eq 

The inverse decay rate is given by 



rv = \r, D + 2 (rjp + r « ) + 2if, + ^r<<> . (59) 
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where and n° q are the equilibrium number densities of heavy neutrinos and lepton 

doublets, respectively, and Ki j2 (z) are modified Bessel functions of the third kind. The 

(x) 

quantities F\ are thermally averaged reaction rates per particle X. They are related by 
= 7i/n^ to the reaction densities 7$ [36]. Our calculations are based on the reduced 
cross sections given in ref. [38]. 

It is very instructive to consider analytical approximations to the various washout 
contributions. Both, the inverse decay rate and the resonance part of F^J (cf. [9]) are 
proportional to Ki(z), 

r # = l r 'D + 2Tg es = ^-L^ (h*h) n M^K^z) . (61) 

The integral in Eq. (56) can be analytically performed. The corresponding washout factor 
can be written in the form 



^>( Zi ) = exp I -2^3) ^ \Y + ** K2 M - Y Zi (KtWHzi) + Hzi)Ki(zi))) j,(62) 

where is the equilibrium mass (2), and L 1>2 (z) are modified Struve functions [49]. 
Rather accurate approximations are, for small and large values of Z\ respectively, 



exp{-Icl3)SH¥ -§*f + W))}, *<1 

exp {"icWt^e- (z? + f , 2 + fg* + fif + },*>!. 



(63) 

The non-resonant contribution of Ni exchange to the washout is proportional to m 2 , 

1 ~~ Z I N,nonres ~^ L N,tJ 

1 iWfm 2 1 . . 

(64) 



7T 3 C(3) V* 1 Z 3 ' 

which yields the washout factor 

f5l/ > f 8 Mim 2 1 I 

- <2, (*>=exp{-^^-}. (65) 

Finally, we have to consider iVi-top scatterings. The rate is dominated by t-channel 
Higgs exchange if the infrared divergence is cut off by a Higgs mass ~ 1 TeV <C T. 
In terms of the reduced cross section one has (cf. [38]), 

^ = 96^(3) / teVx&teW^zy/x) . (66) 
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Figure 8: Washout factor as function of the initial temperature z^ = M\fT\ for different 
values of rh\ and M\ = 10 8 GeV; iVi-top scatterings are neglected. 

For small and large values of z iy respectively, analytic expressions are given by 

ex p{-^S:( ln ( 4 ^)- 1 )}^i< 1 



^ (3) (^i) 



^{-i^ySrv^ 6 "" ( ln (?) (t + *) 



- Zi 



(67) 



Zi > 1 , 



where a u = ^/(47r) and a<f, = Mf/m^. 
The total washout factor 



u(zi) 



N 



(68) 



B—L 



depends exponentially on the parameters rh\ and M{fri 2 {oo^). For not too 

large M x and not too small z\ (cf. Figs. (8)- (10)), oo^ ~ 1 whereas u/ 1 ) reaches a plateau 
for Z{ < 1 at 



u)^\zi) ~ exp 



37r mi 



4C(3)m,y (69) 
At smaller values of Zi, and correspondingly higher temperatures 7], eventually a/ 2 ) de- 
creases rapidly. When 7] reaches M 2 , the mass of N 2 , a new plateau will be reached. The 
larger Mi, the larger the value of z\ where the decrease of sets in. This behaviour is 
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Figure 9: Washout factor as function of the initial temperature Zi = M 1 /T i for different 
values of rh\ and M\ = 10 10 GeV; iVi-top scatterings are neglected. 




Figure 10: Washout factor as function of the initial temperature Zi = Mi/7] for different 
values of rh\ and M 1 = 10 12 GeV; iVVtop scatterings are neglected. 
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Figure 11: Washout factor as function of the initial temperature Zi = M\fT\ for different 
values of rh\ and M\ = 10 14 GeV; Ai-top scatterings are neglected. 

clearly visible in Figs. (8)-(10). At very large Mi, the decrease of oo^ is effective already 
at large values of z\ (cf. Fig. (11)). 

The factor u/ 3 ) is very sensitive to the value of a^, i.e. the choice of the infrared cutoff 
m^. For = 1 TeV, cu^ significantly improves the washout of uj^uj^ 2 \ but it does 
not change the qualitative picture. This is illustrated by Fig. (12) where the cases with 
and without iVi-top scatterings are compared. On the other hand, for ~ Mi, is 
always negligible compared to u/ 1 ). The issue of the correct choice of the infrared cutoff is 
theoretically not yet settled. There is a corresponding, though less important uncertainty 
in the generation of the baryon asymmetry for small values of rh\ [39]. The washout 
factors uj^uj^ shown in Figs. (8)-(ll) can be regarded as conservative upper bounds on 
the full washout factors to = uj^uj^uj^ . 

It is remarkable that the washout of an initial asymmetry at Z{ ~ 1, i.e. 7] ~ Mi, 
becomes very efficient for rh\ > m* ~ 10~ 3 eV. Since the efficiency increases exponentially 
with increasing fhi, already at rhi = 5 x 1CT 3 eV one has uj(z[ = 1) < 10~ 4 . Hence, 
for neutrino masses of order or larger than a/ Am s 2 |, AL = 1 processes are very likely 
to erase any previously generated baryon asymmetry to a level below the asymmetry 
produced by leptogenesis. As shown in [9], for these values of rh\ the final asymmetry is 
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Figure 12: Comparison of the washout factors as function of Zi = M 1 /T i without (full 
line) and with (dashed line) TVi-top scatterings; Mi = 10 10 GeV. 

also independent of the initial Ni abundance. Hence, a complete independence of initial 
conditions is achieved. 

6 Summary 

We have extended our previous work on the minimal version of thermal leptogenesis where 
interactions of N\, the lightest of the heavy Majorana neutrinos, are the dominant source 
of the baryon asymmetry. Based on the seesaw mechanism, we have derived an improved 
upper bound on the CP asymmetry e±, which depends on Mi, the mass of Ni, the light 
neutrino masses m x and m 3 , and the effective neutrino mass m 1 . Given the two mass 
splittings Am^ tm and Ara^, the neutrino masses mi and m 3 can depend on the absolute 
neutrino mass scale m in two ways, corresponding to normal and inverted mass hierarchy, 
respectively. 

From the numerical solution of the Boltzmann equations we have obtained an upper 
bound on all light neutrino masses of 0.12 eV, which holds for normal as well as inverted 
neutrino mass hierarchy. This is about a factor of two below the recent upper bound of 
0.23 eV obtained by MAP [50]. The leptogenesis bound is remarkably stable with respect 
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to changes of r/g MB , Am^ m , the effect of supersymmetry, and theoretical uncertainties 
of 77^ ax . Quasi-degenerate neutrinos are only allowed if the CP asymmetry is strongly 
enhanced by a degeneracy of the heavy Majorana neutrinos. For instance, in order to 
relax the upper bound to 0.4 eV, degeneracies AM 2 i/Mi, AM31/M1 < 10~ 3 are required. 

We have also studied the washout of a large, pre-existing B — L asymmetry. It is very 
interesting that a washout by several orders of magnitude takes place at temperatures 
T close to Mi, if the effective neutrino mass rhi is larger than the equilibrium mass 

~ eV. All memory of the initial conditions is then erased. 

We conclude that for neutrino masses in the range from 10~ 3 eV to 0.1 eV leptoge- 
nesis naturally explains the observed baryon asymmetry, independent of possible other 
pre-existing asymmetries. It is very remarkable that the data on solar and atmospheric 
neutrinos indicate neutrino masses precisely in this range. 
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